Diffractionless propagation through Kapitza stratified media 



Carlo Rizza 1,2 and Alessandro Ciattoni 2 

Dipartimento di Scienza e Alta Tecnologia, Universita dell'Insubria, Via Valleggio 11, 22100 Como, Italy 
Consiglio Nazionale delle Ricerche, CNR-SPIN, 61100 Coppito L'Aquila, Italy 

We show that in the presence of a rapidly modulated dielectric permittivity with a large modula- 
tion depth (Kapitza medium) a novel and robust regime of diffractionless electromagnetic propaga- 
tion occurs. This happens when the mean value to depth ratio of the dielectric profile is comparable 
to the small ratio between the modulation period and the wavelength. We show that the standard 
effective medium theory is inadequate to describe the proposed regime and that, its occurrence is 
not substantially hampered by medium losses. We check the feasibility of the proposed regime by 
means of a large modulation depth metal-dielectric layered medium whose subwavelength imaging 
properties are analytically investigated. 

PACS numbers: 



Subwavelength imaging and diffractionless propaga- 
tion are major photonics subjects that have attracted, 
in the past decade, a renewed research interest triggered 
by the advent of metamaterials [H, 0. In fact, nega- 
tive index materials have been shown to support negative 
refraction of propagating waves and, at the same time, 
magnification of the evanescent waves, thus providing an 
overall superlensing effect. Remarkably, without resort- 
ing to exotic and hardly available magnetic properties, 
it has been realized that a metal slab can provide sub- 
wavelength imaging as well, the effect being supported 
by plasmon excitation at its surfaces [3, 4] . A refinement 
of the metal slab superlens has been proposed in the form 
of metal-dielectric layered composite obtained by alter- 
nating two materials of equal thickness and permittivities 
with real parts of opposite signs [5|,|6|. Basically, the cho- 
sen dielectric-metal symmetrical composition provides, 
within the standard effective medium approach, trans- 
verse and longitudinal permittivities which are zero and 
infinite, respectively, thus effectively enabling the struc- 
ture to support the electrostatic situation where each 
wave is replicated both in amplitude and phase. How- 
ever, due to the high impedance mismatch between such 
structure and vacuum, a different layered metal-dielectric 
composite has been proposed which can be tailored to op- 
erate in the canalization regime where the effective trans- 
verse permittivity is equal to 1 while the longitudinal per- 
mittivity is very large 0413] ■ Such dielectric anisotropic 
response provides the effective medium with a very flat 
equal frequency contour, this allowing the slab to trans- 
port nearly all the spatial harmonics of the impinging 
image and, at the same time, to fulfill the Fabry-Perot 
resonance condition for enhancing transmission. Another 
subdiffraction imaging effect exploits a different kind of 
anisotropic metamaterials whose ordinary and extraor- 
dinary permittivities have different signs (indefinite or 
hyperbolic media) (TT| - [T4| . The related hyperlensing ef- 
fect is based on hyperbolicity of the extraordinary waves 
equifrequency contour whose positive curvature entails 
negative refraction and whose unboundedness intrinsi- 
cally allows propagation of the subwavelength field fea- 
tures. Recently, it has been proposed that a plasma in 



the presence of a strong magnetic field is able to provide 
subdiffraction imaging for the electromagnetic frequency 
well below the plasma frequency [l5j , and also such con- 
figuration exploits the indefinite (anisotropic) response of 
the magnetized plasma whose hyperbolic equifrequency 
contour is very flat at sufficiently low frequencies. 

In this Letter we show that a novel regime of diffrac- 
tionless propagation occurs for Transverse Magnetic 
(TM) waves propagating, in the long wavelength regime, 
through rapidly modulated stratified media with a large 
dielectric modulation depth (Kapitza media). The con- 
sidered electromagnetic situation is conceptually equiva- 
lent to that of the mechanical inverted pendulum whose 
pivot point is subjected to high-frequency vertical oscil- 
lations (Kapitza pendulum) fu| [I?} , since these oscilla- 
tions produce a rapidly varying contribution to the la- 
gragian function with a large modulation depth. The 
concept has been successfully extended and implemented 
in the gene ral context of quantum [Tij and nonlinear 
physics [19| . such as for kink solitons in perturbed sine- 
Gordon models, with applications for long Josephson 
junctions [20| . guiding-center solitons in perturbed Non- 
linear Schroedinger Equations, with applications in opti- 
cal pulse propagation in fiber [2l|, [22j , in quadratic non- 
linearities with applications for spatial quadratic solitons 
[23| and in " nonlinearity management" in both Bose- 
Einstein Condensates dynamics and nonlinear optics (24h - 
I27I ]. In all these cases, even though the system is not 
able to follow the rapid external oscillations, yet these 
are able to affect the average system dynamics by means 
of additional slowly- varying potential contributions (e.g., 
the rapid external oscillations yield the stabilization of 
the inverted state of the mechanical pendulum). In the 
electromagnetic situation we consider in this Letter, the 
rapid and deep dielectric spatial modulation entails the 
diffractionless propagation of the averaged electromag- 
netic field. 

The electromagnetic field amplitudes E = E x (x, z)e x + 
E z (x, z)e z , H = H y (x 1 z)e y associated with monochro- 
matic TM waves satisfy Maxwell equations 

- d x E z + d z E x = iuj[j, Hy, 
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where e = diag[ea;, e y , e z ] is the medium relative dielec- 
tric tensor and the time dependence e~ lut has been as- 
sumed. Here we consider a specific medium periodically 
modulated along the z-axis whose dielectric permittivity 
e x = e z — e admits the Fourier series expansion 
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where e m and (a n + b n /rj) are the Fourier coefficients 
whereas 2nr]/K is the spatial period. Here 77 is a di- 
mensionless parameter we have introduced to explore the 
asymptotic electromagnetic behavior pertaining the limit 
77 — > where both the grating amplitude and its spatial 
frequency are very large (assuming that K w fco = oj/c). 
Note that the coefficients a n account for a contribution 
to the modulation which is not large, and it is essential 
for allowing the model to avoid unrealistic features as a 
permittivity imaginary part rapidly oscillating between 
large negative and positive values. Since electromagnetic 
propagation is here characterized by two very different 
scales (i.e. the radiation wavelength and the permittiv- 
ity modulation period), it is natural to let each electro- 
magnetic field component to separately depend on z and 
Z = z/rj (multiscale technique where Z is the fast coor- 
dinate [28j |) and to represent it as a Taylor expansion up 
to first order in 77, i.e. 



A{x,z,Z)= A i °\x,z) + A < - a) (x,z,Z) 
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where A — E X ,E Z or H y , the superscripts (0) and (1) 
label the order of each term whereas the overline and the 
tilde label the averaged and rapidly varying contributions 
to each order, respectively. Substituting the field com- 
ponents of this form into Eqs.([T]), each equation yields a 
power series in 77 whose various orders are evidently the 
superposition of a slowly varying (independent on Z) and 
a fast (dependent on Z) contribution all of which can in- 
dependently be balanced. From the lowest order (which 



is here 77 1 ) we obtain E : 
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words, the dominant contribution to x-component of 
electric field is slowly varying, whereas the z-component 
of the electric field has no zeroth order. The order 77 of 
the first of Eqs.([T|) yields 



and dzE x ^ — iu>fj,oHy , which in turn, with the help 
of the obtained expression for H y °\ yields Ex — 
(k%/K 2 )J2 n ^ {b (n) /n 2 )e inKZ E x 0) . The slowly varying 
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part of the order 77 of the second of Eqs.(fT]), together 



(1) with the obtained expression for Ex yield 



d z H y = -iuje e e f f E x . 



where 



K 2 



(5) 
(6) 



To summarize we have obtained that the leading term 
of the TM electromagnetic field is slowly varying and 
rigorously Transverse Electro-Magnetic (TEM) , i.e. E ~ 

E x °\x, z)e x , H ~ Hy (x, z)e y , and, in addition, its com- 
ponents satisfy Eqs.(j4|) and ((5|) so that the field rigor- 
ously undergoes diffractionless propagation and experi- 
ences the uniform effective permittivity e e //- 

As for any Effective Medium Theory (EMT), the field 
is here slowly varying since it cannot follow the rapid 
medium dielectric modulation. However, the fact that 
the longitudinal component vanishes is a very unique 
and important feature of the proposed regime which 
has a simple physical explanation arising from the com- 
bination of the first Maxwell Equation and the large 
depth of the dielectric modulation. In fact, a very 
deep and rapidly modulated dielectric response would 
produce volume polarization charges too strong to be 
screened by the electric field, an incompatibility avoided 
by the condition E z ~ 0. For example, assuming e = 
e m + (b/r])cos(Kz/r]), the equation V • D = yields 
e(d x E x + d z E z ) - (K/r) 2 )sin{Kz/r))E z = from which 
it is evident that, since the first term self-consistently 
scales as I/77, E z has to be proportional to 77 for mit- 
igating the divergence of the factor I/77 2 in the second 
term. The diffractionless propagation of the leading con- 
tribution to the electromagnetic field in the considered 
regime is a consequence of its TEM structure since, as 
it is well known, TM waves with vanishing longitudinal 
component do not undergo diffraction (as it is evident 
from Eqs. (TT])). It should be stressed that the super- 
resolution mechanisms reported in literature are based 
on the requirement that e z — > 00 (since, when combined 
with the third of Eqs. (TT]), this condition yields E z — >• 
and hence diffractionless propagation) a condition that 
has to be fulfilled by exploiting the standard EMT and 
which is essentially hampered by medium absorption Q ■ 
On the other hand, in the regime we are investigating 
here, the longitudinal effective permittivity e z is not even 
defined and it is evident from the above presented mul- 
tiscale analysis that the longitudinal field component E z 
vanishes regardless the medium losses (since we have 
not assumed the permittivity to be real). We conclude 
that the occurrence of diffractionless propagation within 
Kapitza stratified media is very robust and substantially 
not hampered by medium losses. Note that, in analogy 
to the potential governing the averaged dynamics of the 
Kapitza pendulum, the effective permittivity of Eq.® 
is the sum of the average permittivity plus a contribu- 
tion arising from the rapidly varying part of the dielectric 
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FIG. 1: (Color online) (a) Layered slab and TM waves scat- 
tering geometry, (b) and (c) real and imaginary parts of the 
dielectric permittivity profile within the unit cell, (d) Com- 
parison between the exact slab transmissivity (solid lines) and 
thoose predicted by the Kapitza (dashed lines) and the stan- 
dard EMTs (dash dotted lines) for a the unit cell dielectric 
distribution of panel (b) and (c) with e m = 0.05 + 0.05i, Se = 
0.025, N = 10 and r\ = 1/60, and wavelength A = 100 \im. 
The semitransparent surface interpolating the standard EMT 
transmissivity profiles is plotted for clarity purposes. 



modulation. Such an effective medium response, together 
with the TEM field structure, clarifies that the Effective 
Medium Theory we are considering here (Kapitza EMT) 
is fundamentally different from the Standard Effective 
Medium theory (Standard EMT) for layered media whose 
main result is that the effective principal permittivities 
are e x = (e) = e m and e z = (e -1 ) -1 . The difference be- 
tween the two regimes is also assured by the fact that the 
standard effective medium results can be simply derived 
by using the above discussed multiscale technique in the 
presence of the dielectric profile of Eq.@ with b n = 0, 
i.e. without the large modulation depth contribution to 
the permittivity. 

In order to check the predictions of the proposed 
Kapitza EMT, we have considered reflection and trans- 
mission of TM plane waves from a slab filled by a Kapitza 
stratified medium, as sketched in Fig. 1(a). The dielec- 
tric modulation is along the 2-axis with period rj\ where 
A = 2n/ko is the vacuum radiation wavelength and 77 
is the above introduced small parameter, whereas the 
slab thickness L is a multiple of the period i]A. The 
unit cell comprises N homogeneous layers of thicknesses 
r)\/N and the dielectric permittivity of the j — ih layer 
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where e m is the mean value of the dielectric permittiv- 
ity and Se is responsible for the (not large) modulation 
of medium absorption. The corresponding profile of the 
slab dielectric permittivity is evidently of the kind of 
Eq.@ (with K = ko) and its real and imaginary parts 
are sketched, within a unit cell, in Figs. 1(b) and 1(c), 
respectively (where e' m = Re(e m ) and ej^ = Im(e m )). 
The check of the Kapitsa EMT has been performed by 
choosing A = 100 fim, e m = 0.05 + 0.05i, Se = 0.025, 
N = 10 and r\ = 1/60, which is a realistic situation 
since -60 < Re(e) < 60 and 0.025 < Im(e) < 0.075 
(thus avoiding the use of active media) and since the 
layers' thickness T)X/N ~ 166 nm is fully feasible. The 
scattering process of TM waves by the considered lay- 
ered medium admits full analytical description (by means 
of the transfer-matrix method) and in Fig. 1(d) we have 
plotted the exact profiles (solid line) of the transmissivity 
defined as T = |E t | 2 /|E;| 2 (see Fig.l(a) for the definition 
of the field amplitudes) as function of the transverse wave 
vector k x — ko sin 9 for various slab thicknesses L. On 
the other hand, Eqs.(|4|) and ([5]) are easily solved to yield 
the transmissivity 



T = \cos{k z L) - iFsm(k z L) 



(8) 



where k z = ko^/e e ff and F = 
[V £ eff cos# + 1/ (y/e e ff cos 0)\ /2. After evaluat- 
ing the Fourier coefficients of the considered dielectric 
profile, Eq.© yields e eff = 0.5339 + 0.05i and in 
Fig. 1(d) we have reported various profiles (dashed lines) 
of the resulting transmissivity of Eq.© pertaining to 
the Kapitza EMT. In addition, the standard EMT 
would describe the slab as an anisotropic medium 
with dielectric permittivities e x = (e) = 0.05 + 0.05i 
and e z = (e" 1 )- 1 = (-7.205 + 7.194i) ■ 10 3 and, in 
Fig. 1(d), we have plotted various profiles (dashed dot 
lines) of the the corresponding transmissivity which 



coincides with Eq.© with k z = k, 



1 - 



and 



F = 



1 / fcpe^ cos 

2 I k z 



r; , . Note that the agreement 

koe x cos & J 

between the exact predictions and those based on the 
Kapitza EMT is remarkable whereas the standard EMT 
completely fails. 

It is well known that medium absorption plays a very 
detrimental role in the observation of the canalization 
regime 8] as well as on all the other mechanisms sup- 
porting diffractionless propagation. In order to investi- 
gate the effect of medium losses on the Kapitza regime 
we are here discussing, we have extended the scatter- 
ing experiment of Fig. 1(a) to encompass even evanescent 
waves and we have evaluated, for a number of slabs char- 
acterized by different absorption efficiencies, the quantity 
T = \Et\1 =L /\Eii\1 =0 as a function of the transverse wave 
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FIG. 2: (Color online) Logarithmic plot of the quantity 
T = |Et|^ =i /|Ei|J = o, as function of the transverse impinging 
momentum k x , generalizing the slab transmissivity of Fig. 1(d) 
to encompass transport of evanescent waves for three slabs of 
thickness L — 1.66A and r\ — 1/60, with different mean values 
and modulation depths of the permittivity imaginary part. 
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FIG. 3: (Color online) (a) Imaging of a subwavelength two- 
peaked shape at A = 100 fim by a Kapitza slab (as reported in 
Fig. 1(a) with unit cell dielectric profile of Figs. 1(b) and 1(c)) 
characterized by e m = 0.05 + 0.005i, 5e = 0.0025, N = 10, r\ = 
1/60 and L — 135.66/im. (b) Diffraction of the subwavelength 
two-peaked shape after L = 135.66 /im vacuum propagation, 
(c) Same as in panel (a) but with r\ — 1/20. 



vector k x (which, for \k x \ < k , coincides with the above 
discussed slab transmissivity and, for \k x \ > ko, entails 
information on the slab efficiency to transport evanescent 
waves). In Fig. 2 we report the exact logarithmic plot of 
T together with the corresponding profiles predicted by 
the Kapitza and the standard EMTs for three different 
slabs of thickness L — 1.66 A which are identical to those 
considered in Fig.l but characterized by different values 
of e m and Se: (a) e m = 0.05 + 0.005i, Se = 0.0025; (b) 
e m = 0.05 + 0.05i, Se = 0.025; (c) e m = 0.05 + 0.5i, 
Se = 0.025. The corresponding dielectric permittivi- 
ties predicted by the Kapitza and the standard EMTs 
are: (a) e eff = 0.5339 + 0.0050i, e x = 0.05 + 0.005z, 
e z = (-14.256 + 1.424i) • 10 3 ; (b) e eff = 0.5339 + 0.05«, 
e x = 0.05 + OMi, e z = (-7.205 + 7.19U) ■ 10 3 ; (c) e eff = 
0.5339+0.05i, e x = 0.05+0.5i, e z = (-0.154+1.425i)-10 3 . 
From Fig. 2 we note that in all the three considered sit- 
uations the Kapitza EMT predictions largely agree with 
the exact phenomenology whereas the standard EMT is 
inadequate and, most importantly, that this happens al- 
most regardless the medium absorption. Considering that 
in panel (c) of Fig. 2 the mean value and the modulation 
depth of the permittivity imaginary part are one hundred 
times greater than those of panel (a), we restate that the 
Kapitza regime is very robust against medium losses, a 
very distinct and important feature which does not char- 
acterize the other diffractionless propagation regime pre- 
sented in literature. 

From Fig. 2, it should also be stressed that the agree- 
ment between the Kapitza EMT and the exact phe- 
nomenology is, in the considered situation, very good and 
robust even for evanescent waves ( | fc^ | > ko), an essential 
requirement for achieving super-resolution. Even though 
we have mainly focused on the novel diffractionless prop- 



agation regime, it is worth checking the subwavelength 
imaging features of a Kapitza layered slab. Specifically 
we have considered a slab whose unit cell has the di- 
electric distribution of Eq.dTJ) with e m = 0.05 -I- 0.005i, 
Se = 0.0025, N = 10 and r) = 1/60 (thus coinciding 
with the case of Fig. 2(a)). The radiation wavelength is 
A = 100 iim and the unit cell thickness is A = t]\ = 
1.67 /im. We have considered a slab with M — 82 unit 
cells since its resulting width L = MA = 135.66 fim is 
very close to the width Lp = A/^/Rc(e e y^) = 136.85 /im 
for which the Kapitza homogeneous slab governed by 
Eqs.g]) and © with e eff = 0.5339 + 0.005i (from Eq.©) 
has its second Fabry-Perot resonance. Using the transfer- 
matrix method, we have evaluated the slab point-spread- 
function PSF(x) according to which the transmitted 

magnetic field profile Hy T ^ (x, L) can be expressed as 
the convolution of PSF{x) with the incident magnetic 



field H^\x,0), 



or H^\ Xl L) = PSF(x) * H^'{x 7 0) 



30]. In Fig. 3(a) we plot the comparison between the 
slab output and input images for an impinging magnetic 
field whose profile consists of two Gaussians of equal 
width a = A/20 whose separation distance is 4er, i.e. 



HW(x,0) = Ho 



-(x-2<j) 2 /<t 2 



Apart 



from a slight decrease of the peaks height due to medium 
absorption, the outcoming image is clearly a very ac- 
curate copy of the incoming one. To appreciate such 
imaging result, we have reported in Fig. 3(b) the effect of 
vacuum propagation on the above incoming field, and the 
complete deterioration of the considered subwavelength 
image is evident. For completeness, in Fig. 3(c) we have 
plotted the analogous fields comparison for a slab identi- 
cal to that of panel (c) but with 77 = 1/20 whose subwave- 
length imaging efficiency is evident even in the present 
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of a slight image deterioration resulting from the small 
sidelobes. 

In conclusion, we have proposed and discussed a novel 
diffractionless propagation regime which is not funda- 
mentally affected by medium losses and it is very ro- 
bust against the tailoring of the medium dielectric prop- 
erties. Diffraction suppression is here achieved through 
propagation in a medium whose dielectric permittivity 
is rapidly modulated with a large modulation depth so 
that the physical underpinnings of the proposed setup 
arc in many respects an elaboration of those that in me- 
chanics are associated with the stabilization of the in- 
verted pendulum (the Kapitza pendulum from which the 
name Kapitza EMT for our approach) . In addition to the 



broad interest of our results (obtained from the physical 
ubiquitous situation where a rapidly varying and strong 
external stimulus is able to affect the slow system dy- 
namics) we stress that our method is amenable to many 
relevant elaborations, as for example, the extension to 
the case where also the medium magnetic permeability 
is rapidly varying with a large modulation dept. We sug- 
gest that such a situation, surpassing the here consid- 
ered TM electromagnetic case, could in principle be able 
to support diffractionless propagation of the full electro- 
magnetic field. 
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